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We analyze the singlet P-wave charmonium production at e + e~ colliders within the framework 
of unconstrained supersymmetry. We show that the CP-violating transitions, dominated by the 
gluino exchange, are typically four orders of magnitude larger than the CP-conserving ones, and 
former is generated by the electric dipole moment of the charm quark. Our results can be directly 
tested via the charmonium searches at the CLEO-c experiment. 
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§ ' I. INTRODUCTION 

CN ■ 

^ ' A permanent electric dipole moment (EDM) of an elementary particle is a clear signature for the CP violation. 
In the Standard Model (SM) of electroweak interactions CP-violation originates from the phase in the CKM matrix; 
but in this model the single fermion's EDM vanishes up to two-loop order, and also three-loop contributions partially 
cancel. Therefore, the SM contribution is neither relevant phenomenologically nor numerically important. However, 
^1 ■ the Supersymmetric (SUSY) extension of SM contains many new sources of CP violation generated by soft SUSY 
breaking terms. Actually such a SUSY Model has many phases, but not all of them are physical. After using some 

■ symmetries such as R-symmetry and the Peccie-Quinn symmetry, only two of them remain physical: the trilinear 
^ ' coupling A and the fi parameter. 

The upper bounds for the neutron and electron EDMs imply that CP violation phases are duly small if the exchanged 
SUSY particles have masses close to their current experimental lower limits. There are several suggestions to find 
the electron and neutron EDMs below the experimental upper limits when CP phases ~0(1). Two possibilities have 
been commonly discussed in the literature which include (i ) finding appropriate parameter domain where different 
\ contributions cancel or (ii)making the first two generations of scalar fermions heavy enough but keeping the soft 

■ masses of the third generation below the TeV scale. 
As the analysis of Ref. [1] shows clearly, the electric dipole moments (EDM) of the heavy quarks play an important 

role in the direct production of singlet P-wave mesons at e + e~ colliders. Indeed, it is possible to observe x Pi 
bottomonium in B factories like BaBaR, KEK or BEPC with sufficiently high statistics provided that the b quark 
possesses a large enough EDM. Since the single fermion EDMs are exceedingly small in the standard electroweak 
theory [2] the SUSY contribution remains as the only viable option. Indeed, in SUSY the EDMs exist already at the 
one-loop level and typically exceed the existing 1.5a upper bounds. 

However, not only the bottom quark but also the charm quark EDM offers an important arena in searching for 
the singlet P-wave charmonia. Indeed, in near future with increasing data at CLEO-c experiment [3], it might be 
possible to produce the 1 Pi charmonium state, h c , directly. Therefore, a direct estimate of the charm quark EDM V c 
in SUSY will help in predicting how large the effect could be at CLEO-c. 

In our analysis we will concentrate on unconstrained low-energy SUSY model as in Ref. [1] in that portion of the 
parameter space where large SUSY contributions to electron and neutron EDMs are naturally suppressed. Such a 
region of the SUSY parameter space is characterized by a small phase for the /x parameter: — 7r/(5tan/3) < ip^ < 
7r/(5 tan /?) [4]. 

Assuming universality of the gaugino masses at some ultra high scale, it is known that the only physical soft phases 
are those of the A parameters, ifA { (i = u, c, ■ ■ ■ b) and the fi parameter. Therefore, in what follows we will scan the 
SUSY parameter space by varying (i) all soft masses from m t to 1 TeV, (ii) the phases of the A parameters from 
to 7T, and finally, (Hi) the phase of the \i parameter in bounds noted above suggested by the cancellation mechanism. 

In the second section, we calculate the chargino and gluino contributions to the c-quark EDM in effective SUSY at 
one-loop level and show that the gluino contribution dominates. 

Section 3 is devoted to the analysis of our numerical calculations. 

In section 4, we discuss the possible signatures of the c-quark EDM in the e + e~ annihilation. After examining 
the formation of 1 P\ charmonium state in e + e~ channel, we show that the CP-violating part, generated by the c- 
quark EDM dominates over the CP-conserving ones. In the case of observing the CP-odd resonance in the data 
accumulating at CLEO-c experiment, this would be a direct evidence for the charm quark's EDM and also for the 
CP- violation in SUSY. 
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2. THE CHARM QUARK EDM in SUSY 



The EDM of a spin 1/2 particle is defined as follows 



(1) 



which obviously appears as a quantum loop-effects. In the SM, Vf vanishes up to two- loop order [5], whereas in SUSY 
it exists already at one-loop level. 

As shown in Fig. 1 the c EDM is generated by the gluino, chargino and neutralino diagrams at one loop level. The 
general expressions for the quark EDMs in SUSY can be found in [6] . At large tan f3 the individual contributions to 
the c quark EDM take the following compact forms. The gluino contribution is given by 
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where the loop function B(x) can be found in [6]. Here M~ ± are the masses of the charm squark. It is clear that at 
large tan /3 the (dominant) gluino-scharm contribution to c-quark EDM is almost independent of the phase of the \i 
parameter. 

The chargino contribution, which proceeds via the exchange of charged gauginos and Higgsinos together with the 
scalar s-quark does have a direct dependence on the phase of the /i parameter via chargino mass matrix [1]. The 
exact expression reads as 
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where the first line results from the direct computation, and depends on the vertex factors Vj± and the loop function 
F± both defined in Ref. [1]. The second line, however, is the dipole coefficient arising in c — > uj decay as was pointed 

out in Ref. [1] for b — ► s"f. (This relation to rare radiative decays is particularly important as Cy (M\y) is particularly 
sensitive to large tan j3 effects [7] which introduces large enhancements compared to the SM.) A precise determination 
of the branching ratio and the CP-asymmetry in c — > u7(or the hadronic mode D — > 777 ) [8] decay will help in fixing 
the chargino contribution to the c-quark EDM. Presently, the experimental data is not precise enough to bound this 
contribution so that we will analyze the result of the direct calculation in (3) in the numerical studies below. 

The neutralino contribution is at least 20 times smaller numerically than the gluino contribution. So we do not 
include it in our numerical analysis below. 

3. NUMERICAL ANALYSIS 

We now discuss numerically the gluino and chargino contributions to the EDM of the c-quark. As mentioned before, 
in our analysis we take the SUSY parameters in the following ranges: 

• Soft masses: m| < M? M? \A C . S \ 2 < (1 TeV) 2 



Gaugino Masses: m t < M 2 < 1 TeV with m~ g = fjM 2 
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• Phase of the ft parameter: < p„ < 5^ 

• Phase of the A parameters: < fA CiS < ^ 

• tan/3: 10 < tan/3 < 50 

where the modulus of fi is determined via the relation 

1 



\A C - M cot 0\ = ^ [(Ml - M?f - rM*, (cos 2/3) 2 j (4) 

with r = {gl (5/3)flf)/(2 ff |). 

In Fig. 2 we plot the gluino contribution, D*?, in units of 10 26 e.cm, as a function of the A c phase (Fig. 2(a)), \i 
parameter phase (Fig. 2(c)) and tan/3 (Fig. 2(b)). We have to emphasize that here we use the exact formula instead 
of Eq.(2). In all these plots, we scan the other parameters in their allowed ranges. It is clear from Fig. 2(a) that for 
most of the tp A space (i.e. from to it) the value of the EDM is around 10 -21 e.cm, showing a slight maximum at ir/2. 
At (fA—0, the 10 2 times smaller value of T> 9 C comes from the ip^ phase in the large tan/3 limit. As we have already 
mentioned, Fig. 2(b) shows the almost independent behavior of the EDM on ip^ phase for large tan/3. In Fig. 2(c) we 
see a slight dependence on tan/3 , as we expected for the c-quark contrary to the case of the b-quark [1]. 

We conclude that for a reasonable portion of the parameter space the c-quark EDM gets a gluino contribution to 
the order of 10~ 21 e.cm. 

Fig.3(a-c) show the dependence of the chargino contribution on the parameters <pa, and tan/3, respectively. All 
the plots give a very small contribution, smaller than 10~ 24 e. cm which is at least two orders of magnitude smaller 
than the gluino contribution. Although its numerical value is far below of the gluino contribution, it would be worth 
mentioning the rather slight increase of the chargino contribution with ip^ and tan/3. In these plots we used exact 

expression for V} which depends on the phase of \i parameter and tan/3 in a very complicated way. In order to 

explain the origin of these behaviours, one needs an analytical dependence of J2t=i Y^j=i ^ ^± factor on ip^ and 

tan/3 which could be achieved only by using some approximations. 



4. CHARM QUARK EDM and 1 PiCHARMONIUM 

Recently the E760 collaboration [9] at Fermilab have announced the first possible observation of the : Pi CP-odd 
charmonium state in pp annihilations and the mass value of h c (lP) is given (3525,20±0.15 ± 0.20) MeV. However 
the production of the 1 Pi charmonium state in the e + e~ annihilation is very interesting from the point of view of 
the experimental evidence of the c-quark EDM, since the coupling of photon to the h c (l 1 Pi (l -1 )) charmonium is 
identical to the effective Lagrangian (1). The quantum numbers of this CP-odd resonance coincide with those of the 
current density [10] 

J a (cc\ 1 P 1 ) = c(x) a -f 5 c(x) (5) 

In the e + e~ annihilation the ^istate can be produced via the jZ and ZZ box diagrams in the framework of the 
SM. After using the same consideration in Ref. [1] one can reach the following effective CP invariant Hamiltonian 

H SM = ■^^G F m e m c <BJ a (cc\ 1 P 1 ) . ^(e+e-^Pi) (6) 

where the current is defined in Eq.(6), and the box function 25 comes from standard loop integrals [11] and it behaves 
as 



*~T7^M— ) (7) 

In minimal SUSY with two Higgs doublets, in addition to jZ box diagram there is another CP-odd Higgs scalar, 
A which also contributes to the formation of 1 Pi resonance in e + e~ annihilation. Replacing the Z boson by A in 
Fig 4, the SUSY contributions to the CP-conserving effective Hamiltonian can be written 
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Hsusy = H S m[M z < — > M A o] 



(8) 



which is typically smaller than the SM amplitude for M^o > Mz- In this sense size of the CP-conserving transitions 
is fixed by the SM not by the SUSY contribution. This is an important difference between the charm and bottom 
Pi production in e + e~ collisions. 
In addition to the CP-conserving decay modes, the 1 P\ state can also be produced via the c-EDM in e + e~ channel 
which is a CP-violating mode. Grey blob in Fig. 4(b) represents the effective Lagrangian (1). In this CP-violating 
mode, e + e~ system is in the most probable state 3 Si. Therefore the effective Hamiltonian violating CP from the Fig. 
4(b) is 



Hsusy = -vjt(~ ) J a(cc\ Pi) . (e + (x) j a e~(x)) 



(9) 



which is a pure SUSY effect. The comparison between the sizes of this CP-violating transition with the CP-conserving 
one is a comparison between the SUSY and SM contributions. In fact, the CP-violating production amplitude 
dominates the CP-conserving one provided that the EDM of charm quark exceeds the critical value 
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which is one order of magnitude smaller than the one for the b quark EDM. If T> c is large enough compared to 
T>c rlt then the CP-violating transition Eq. (9) can be suited to generate the h c meson. In the preceding section we 
performed a scanning of the SUSY parameter space within the bounds mentioned before, and found that V c is well 
above the critical value Eq.(10). 

Now we discuss our estimates with possible h c signatures at charm factories. The bottomonia and charmonia 
production in lepton and hadron colliders have been the primary step for experimental investigation of B and D 



mesons. Indeed, several experiments like the B factories (BABAR, CLEO and KEK-B 



colliders running at 



T(4S), the D factories (CLEO-C and BES e e + colliders running at J/ip, ip', ip" and ^(4140) resonances) as well as 
pp (FNAL E789 at y/s = 800 GcV, and CERN WA102 at ^fs = 450 GcV) and pSi (FNAL E771 at VI = 800 GeV) 
colliders form an arena where different bottomonia and charmonia can be produced in intermediate and final states. 

Among all bb and cc resonances we are particularly interested in the singlet P-wave states, hb{nP) and h c (nP), as 
these are CP-odd states their production and decays are highly sensitive to sources of CP violation in the underlying 
theory Presently, there is no experimental evidence for such resonances except for the discovery of h c (lP) at FNAL 
E760 experiment [9]. However, as the recent work [1] suggests, the production of the singlet P-wave mesons (of 
bb type) at e + e~ colliders are particularly interesting in that there are spectacular enhancements in CP-violating 
transitions if the underlying theory admits large enough electric dipole moments (EDM) for heavy quarks. In this 
sense, production of such resonances is particularly sensitive to the CP violation sources beyond the standard CKM 
picture as the single fermion EDMs in the standard model (SM) are highly suppressed (See [2] and references therein). 
From the experimental point of view, e + e~ colliders running at yfs ~ 2mb (\/s ~ 2m c ) form the basic environment 
where hb(nP) (h c (nP)) resonances can be directly produced. Therefore, the B (e.g. BABAR) and D (e.g. CLEO-C 
[3] and [13]) factories provide the requisite experimental opportunities. Unlike the B factories where the e + e~ center- 
of-mass energy is fixed to T(4S) resonance, in D factories, especially CLEO-C [3], has a variable center-of-mass 
energy enabling one to wander energies around \fs ~ 2m c where h c (\P) can be directly produced. 

Considering CLEO-C for definiteness and its variable center-of-mass energy, it is possible to estimate to what 
extent it can observe h c meson. Let us suppose that CLEO-C observed resonance at or in close vicinity of the E760 
value of Mh then the cross section for this event will be 
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where Rp and R$ are the wavefunctions of h c and 3 Si states, respectively. Clearly, when V c /e <~ M^ 1 the two 
cross sections will be of similar size making h c highly observable. However, except for FNAL E760, no experiment 
has been able to detect such a CP-odd meson so far. Therefore, its production cross section must be rather small 



compared to a(e + e 



Si). This rareness of h c production is directly tied up to the smallness of its production 
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cross section. Indeed, for the face value of T> c /e ~ 10~ 21 cm as follows from the dominant gluino production, one 
expects cr(e + e~ — > /i c ) <~ 10~ 12 cr(e + e~ Si). If CLEO-c observes this resonance with a much larger rate then one 
will need other sources of CP violation to enhance the charm quark EDM. On the other hand, if the observed cross 
section turns out to be this size then it will clearly be an indication for SUSY CP violation. 



5. LIGHT GLUINO CASE 



It is useful to look for regions of the SUSY parameter space where h c production cross section in other words 
the EDM of charm quark is substantially enhanced. This would be useful, at least, for excluding certain portions 
of the SUSY parameter space depending on the size of the signal in near-future experiments like CLEO-c. It is 
straightforward to observe that the charm quark EDM, especially the gluino contribution, is grossly enhanced in 
regions of the SUSY parameter space where the gluino and lighter scharm are nearly degenerate and close to the 
hadronic scale (See [14] for analogous discussion of the scalar bottom quarks). To prevent any conflict with the LEP 
results, we must suppress Z boson couplings to scharms [15] hence we take 



sin #j = y- sin#vv, (12) 
V *j 

and 3> Mj 2 . Here 9c is the mixing angle of scalar charm quarks, and its is fixed by requiring that Z boson does 
not couple to light scharm, c 2 . On the other hand, the two scalar charms must be well splitted in mass to suppress 
e+e~ — > Z* — ► c 2 c 2 and c 2 ci signals at LEP. Under these conditions, the gluino contribution to the charm quark 
EDM becomes, 

V A 9 ~ a s( m g) si rnp Ac ^ 



e J 27n trig 

when m,g ~ Mj 2 . Therefore, when rrig value lies between and 2mt the charm quark EDM is around 
10~ 17 cm. smLpA c which is nine orders of magnitude larger than the critical value computed in Sec. 4. It is clear that, 
with such an enhancement the h c production cross section obeys roughly a(e + e~ — > h c ) w 10~ 4 er(e+e~ — > 3 Si) which 
must be in the experimentally detectable range. With a heavy SUSY spectrum, as investigated in previous sections, 
the h c production cross section is rather small, and unless the detector has a good energy resolution it is difficult to 
observe an h c resonance in s channel. However, a SUSY spectrum with light gluinos and scalar charm quark (with no 
conflict with the existing experiments) it is possible to enhance the h c production cross section significantly making 
it easier to search for such resonances at, for instance, the CLEO experiment. 

For making our results more transparent to experimental investigation, it is useful to discuss the quantity 



R = La(e+e- -> h c ) ~ 27 x N 



M h V c 



J/4> 

where Nj/^ w 10 9 is the expected number of events at CLEO-c. For a heavy sparticle spectrum, R <~ 10 -2 whereas 
for a scenario with light gluinos and scharms one has R ~ 10 5 which is well inside the observable range provided that 
10 9 J/ip mesons are produced. 



6. DISCUSSIONS AND CONCLUSION 



We have discussed the production of singlet P-wave charmonia at e + e~ colliders. This process is determined 
mainly by the CP-violating transitions shown in Fig. 4 (b). It is clear that the EDM of the charm quark is the 
basic machinery to generate h c , and it is only via SUSY with explicit CP violation that the CP-violating transition 
dominates the CP-conserving channel. 

The channel we discuss here, which can be directly tested at the CESR-c collider, suggests a direct access to the 
SUSY phases. The amount of CP violation provided by SUSY is large, and hard to realize in other extensions e.g 
two-doublet models. Therefore, any evidence for h c at the CESR-c experiment will form a direct evidence for SUSY 
in general, and SUSY CP violation in particular. 
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FIG. 1. The Feynman one-loop diagrams of the SUSY contributions to the EDM of c-quark. 
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FIG. 2. The gluino contribution to c-quark EDM in units of 10 26 e.cm as a function of (a) trilinear couplings A c phase, 
fj, parameter phase and (c) tan/3. 
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FIG. 3. The chargino contribution to c-quark EDM in units of 10 26 e.cm as a function of (a) trilinear couplings A s phase, 
(b) /i parameter phase and (c) tan/3. 
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FIG. 4. The Feynman diagrams of the 1 Pi charmonium resonance in e + e scattering: (a) The CP-conserving decay mode 
and (b) is the CP-violating decay mode where the grey blob stands for the insertion of the effective Lagrangian (1). 
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